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I. INTRODUCTION 

A host of cosmological measurements indicate that the 
universe is undergoing a phase of accelerated expansion. 
This has been generally attributed to a significant com- 
ponent of smooth energy with a large negative pressure, 
referred to as dark energy (DE) and characterized by an 
equation of state parameter w = ^. Current measure- 
ments indicate that about 70% of the density of the uni- 
verse today is comprised of this dark energy. Candidates 
for DE include a cosmological constant A, and a slowly 
evolving dynamical scalar field such as quintessence 
In quintessence models, the cosmic acceleration is driven 
by a scalar field (fi slowly evolving in some potential V{4>). 
In this scenario, the parameters of the potential ¥{(/)) de- 
termine the properties of the dark energy. 

In general all DE models have serious unresolved the- 
oretical problems, and one can make the case in different 
ways as to which types, if any (i.e. A or quintessence 
DE), are best motivated 0, l^]- This paper is motivated 
by the fact that scalar field quintessence is definitely part 
of the theoretical discussion, and thus it should also be 
part of the process whereby we evaluate future dark en- 
ergy experiments. This paper is 5th in a series of papers 
motivated in this way [1, S Q • The IPL model we con- 
sider here is one of the more popular quintessence models. 
One of its attractive features is its "tracking" behavior 
that make its predictions independent of the initial con- 
ditions for 0, assuming that (f) starts out in the (rather 
broad) basin of attraction for tracking. Also, the behav- 
ior of the equation of state in the IPL model tends to 
be quite different than for the models considered in our 
previous work (see Q for a unified discussion), so this 
makes it an interesting complement to our other work. 

Recently, the Dark Energy Task Force (DETF) pro- 
duced a report that considered the impact of various pro- 
jected data sets (referred to as "data models" and repre- 
senting future DE observations) on cosmological param- 



eters in a standard ACDM cosmological model using the 
"ifo — Wa" parameterization of the dark energy equation 
of state Q, w(a) = wo + Wa{l — a), where the scale factor 
a — 1 today . They assessed the impact of a given data 
set using a "Figure of Merit" (FOM) , defined as the in- 
verse of the area inside the 95% confidence contour in the 
wq — Wa plane for a fiducial ACDM model. However, as 
has been pointed out by a number of authors (e.g., [lO|), 
the two- parameter wq — Wa phenonienological model is 
not motivated by an actual physical model of dark en- 
ergy and exhibits very different behavior compared with 
popular dark energy models. Our work (represented by 
this and our companion papers 0, [E @) supplements 
the work of the DETF by assessing the capability of fu- 
ture experiments to constrain DE by using an equation 
of state parameterization that is motivated by a physical 
model of DE - the well-known inverse power law (IPL) 
or "Ratra-Peebles" (RP) quintessence model. This po- 
tential has its own motivations, and is also included here 
because it generates a family of functions 'w{a) that are 
quite different than those considered in our other work. 

This paper is organized as follows. In Section[TT]we de- 
scribe the features of the IPL quintessence model and its 
tracking properties. While most of the focus of this pa- 
per is on the tracking behavior of the IPL model, we also 
briefly discuss the non-tracking transient and "thawing" 
behaviors of this model. In Section IIIII we describe how 
we parametrize the IPL model for our MCMC analysis. 
In Section flVBI we present our MCMC analysis and re- 
sults using data forecast by DETF to constrain the IPL 
quintessence model around a fiducial ACDM model. In 
Section flV CI we give our MCMC analysis for simulated 
data generated from a fiducial IPL model. This allows 
us to further ascertain how sensitive future observations 
may be to deviations from a cosmological constant and 
to assess to what extent we can exclude the A model if 
IPL quintessence occurs in nature. In Section IIVDI we 
briefiy discuss our MCMC analysis of non-tracking re- 
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gions of parameter space. Finally, we discuss our results 
and present our conclusions in Section fVl 



II. TRACKING QUINTESSENCE 

For a homogeneous scalar field in an FRW universe, 
the evolution of the scalar field, given by its equation of 
motion, is described by the Klein-Gordon equation 



(1) 



where the Hubble parameter H is given by the Friedmann 
equation (with and spatial curvature also taken into 
account here) 
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where a is the scale factor, Mp = 8ttG^^ is the reduced 
Planck mass, Pr{a) is the radiation background energy 
density, pm{a) is the matter background energy density, 
Pif,{a) is the scalar field energy density, and k is the cur- 
vature constant. The energy density and pressure of the 
scalar field are 
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where the dots denote derivatives with respect to time. 
Eq. (HJ-Eq. ([4]) enable us to solve for the background 
evolution in a quintessence cosmology, once the poten- 
tial V{(l)) and energy densities of the different compo- 
nents, prm Pn etc., have been assigned. If the scalar field 
rolls slowly enough such that the kinetic energy density 
is much less than the potential energy density, i.e., the 
slow-roll limit, 0^ << V{(j)), then the pressure of the 
scalar field will become negative and the field energy will 
approximate the effect of a cosmological constant. This 
indicates that a flat potential V{(j)) is required to give 
rise to accelerated expansion This slow-roll limit 

corresponds to — ~1 and = const. It also follows 
that the equation of state of quintessence is bounded in 
the range — 1 < < 1 and is usually non-constant. In 
these models, the dark energy behaves as a perfect fluid 
in which the equation of state 
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changes with time and is typically negative when V{(j)) 
is sufficiently dominant, as expected during the recent 
epoch of accelerated expansion. We can see from Eq. ([5]) 
that (/> = corresponds to the limit in which the scalar 
field is a cosmological constant with w\ = — 1. 



A. Tracking Solutions and behaviors 

It has been demonstrated [l^, [l^ that a subclass of 
quintessence potentials, including the IPL potential, have 
several desirable properties. These include the fact that 
the equation of motion of these quintessence models have 
attractor-like solutions in the space of trajectories of (f> 
(called "tracking" solutions) . A broad set of initial con- 
ditions and in the early universe (referred to as a 
" basin of attraction" ) evolve toward a common attrac- 
tor solution giving the same late time evolution of 0, 
and thus allowing the scalar field to induce the present 
phase of accelerated cosmic expansion starting from a 
large range of initial conditions. The tracking solutions 
are characterized by an almost constant w^, constrained 
by — 1 < < wb, where is the equation of state of 
the dominating background fluid component. The track- 
ing behavior allows the value of the accelerating mat- 
ter density today to be determined by parameters in the 
quintessence potential, largely independent of the scalar 
field initial conditions We note, however, that al- 

though this behavior may help to explain why the dark 
energy has come to dominate in recent times rather than 
some earlier epoch, it does not solve the "cosmological 
constant problem", especially as it relates to the zero 
point energy of the quantum vacuum. 

In uM. a function 
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(where the primes denote derivatives with respect to (jy) 
was defined for determining whether a particular poten- 
tial admits tracker solutions. It was shown that tracking 
behavior occurs when either of the following two con- 
ditions are met: (a) F > | , wj, < wg , F w const, 



(and thus 
^(1 
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decreases as V decreases); or (b) F < I 

Wb) > Wcf, >wb,F« const, (and thus 

strictly increasing as V decreases). The only constraint 
on the initial energy density in the tracker is that it be 
less than or equal to pb,i, the initial energy density of 
the background fluid component (matter or radiation), 
and greater than Pm.o, the current matter energy den- 
sity. This condition is necessary in order for (j) to converge 
to the tracker solution before the present time [l^ ■ 
On the other hand, solutions of the Klein Gordon equa- 
tion do not converge to tracker solutions for potentials in 



strictly increases as V decreases 
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. Note that 
V 



which Wif, < Wb and 

(F < 1), or, equivalently, when F < 1- 6+2tuB 
^ gives the slope of the potential. The quantity is 

also known as a "slow- roll parameter" (e.g., [9|) which 
relates to how fast the field moves in the potential for 
so-called "slow roll" solutions. One upshot of the above 
analysis is that one can see that potentials (such as IPL) 
tend not to have tracking solutions when and where they 
are fiat (that is where V" = V' = 0). 
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B. The Inverse Power Law Potential 

One of the earliest proposed, simplest, and most widely 
investigated of the scalar field quintessence models is the 
pure inverse power law (IPL) model, originally intro- 
duced by Ratra and Peebles This model was orig- 
inally put forward to mimic a time-varying cosmological 
constant undergoing dissipationless decay and is moti- 
vated by supersymmetric QCD (see [l^ and references 
therein)^ More recently, this potential has been reana- 
lyzed (|12, llSj) in the context of a scalar field potential 
driving the current epoch of cosmic acceleration. 

The IPL scalar field potential is self-interacting, mini- 
mally coupled to gravity, and given by 



(7) 



Values of Vb of order the critical density pc = ZH^Mp 
and a — 0{\) yield cosmological solutions in which the 
scalar field can account for the observed cosmic accelera- 
tion today (and typically has current values = 0{Mp)). 
Furthermore, a large range of cosmologically realistic so- 
lutions exhibit "tracking" behavior whereby, after some 
initial transient period, many different solutions lock on 
to the same attractor solution. This causes the initial 
conditions for 4> to be irrelevant for predicting observ- 
able cosmological features and removes the need for tun- 
ing of initial conditions seen in many other quintessence 
models. 

It has been shown that the following relation is main- 
tained on the attractor solutions [l|, [12, [I3| : 
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The second derivative of the potential gives the scalar 
field mass which today is given by — V^"((/)o) ~ 
The tiny value of this mass (m^ ^ IQ'^^eV) is due to 
the requirements that V{(l)) slowly varies with the field 
value and that the current value of V{(t>) be consistent 
with observations Tz^. When the scalar field potential is 
about to dominate we have using Friedman's equation, 
~ -jpr- Then, if and a are of order unity, Eq. ([8|) 
indicates that the value of the quintessence field at the 
present time is of order of the Planck mass[l8|. 

The power law index a > determines the shape of 
the potential as well as the value of today. The slope 
and curvature of the IPL potential are given by 
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and 
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We can see that smaller a's lead to a more flat potential 
which will in turn lead to more slowly evolving behavior 



for (j) (and thus values of closer to —1). Larger values 
of a lead to a steeper potential slope, causing more evo- 
lution for (j) and its energy density and also values of 
larger than —1. 

Smaller values of (f>i as well as larger values of a lead 
to a steeper initial potential slope and larger values of 
V{(t>i). This means that the scalar field will start rolling 
from higher up on the potential and will roll faster, even 
for cases where the dark energy is initially dominant and 
a is correspondingly large, leading to greater evolution 
of the dark energy density. The quantities Vq and a are 
the two free parameters in the potential. In some su- 
persymmetric QCD realizations of the IPL model [l^ , a 
is also related to the number of flavors and colors, and 
can take on a continuous range of values a > [l9| . 
For a ^ oo (but with still subdominant), the scalar 
field energy density scales like that of the dominant back- 
ground. Potentials of this type also possess the following 
phenomenological property: they yield values which 
automatically decrease to negative values at the begin- 
ning of matter domination [20|. Given that the energy 
density of each component evolves as 



P i oc a 



-3(l+-!«i) 



(11) 



(with i standing for the radiation, matter, or scalar field 
component), quintessence will eventually come to dom- 
inate the universe even if it begins as a subdominant 
constituent. 

The IPL potential is one of a large class of quintessence 
models with what has been referred to as "runaway scalar 
fields" [13, whose tracker solutions begin from some 
initial and (pi and share some of the following gen- 
eral features: The field rapidly converges to a point on 
the potential where V" ~ H^, where the Hubble pa- 
rameter H is determined by pm and pr- As the universe 
expands and H decreases, (j) moves down the potential 
so as to maintain the condition V" « H^. The universe 
enters a tracking phase where catches up to the back- 
ground density ps when decreases to of order 
and so (po ^ Mp [HI, [l3]- Thus, the distinctive feature 
of these tracker fields is that the evolution of the scalar 
fields is controlled by pm and pr rather than evolving 
independently according to its own potential. This con- 
trolled evolution continues until (j) finally surpasses the 
point where critical damping via Hubble expansion oc- 
curs. Then the field's own potential energy is sufficient 
to freeze the field and cause p^ to eventually overtake 
Pm and driving the universe into a phase of cosmic 
acceleration. 

Figured] illustrates how the shape of the IPL potential 
is changed by selecting four different a values for a fixed 
Vq. The value of 4>i determines where on the potential 
the scalar field starts to evolve. The present field value 
00, of order of the Planck mass Mp, is reached from 
a broad range of initial conditions 0/ and 0/, with the 
only important condition being that << Mp [2l| . 
as consistent with the discussion concerning tracking in 
Section FlI Al and the more detailed discussion and criteria 
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regarding attractor solutions given in ll3| . The lower 
panel of Fig. [1] shows the corresponding evolution of 
the equation of state. For fixed values of 0/ and Vq, we 
see that larger values of a correspond to w curves with 
larger amplitudes and which have larger values today, i.e., 
deviate more from a cosmological constant {w = —1) at 
the present time. As a — > 0, the equation of state more 
and more mimics the behavior of A at late times with 
w — » — 1. The IPL model has been categorized by [l^l as 
a "cooling" or "freezing" model in which w > —1 initially 
but with w then decreasing towards — 1 as the scalar field 
rolls down the potential. 
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FIG. 1: IPL potentials (top panel) and w{z) evolution (lower 
panel) for different a values (dashed-dotted: a — 0.05, 
dashed: a = 0.01 and solid: a — 0.1). For all curves Vb ~ 0.38 
and — 10"^". Smaller values of a lead to flatter potentials 
and smaller V{4>). 

For cases in which radiation or matter are dominant 
and the contribution of to the expansion of the uni- 
verse is neglected, the Klein-Gordon equation gives ex- 
act tracking solutions for the evolution of <j) for the IPL 
model, as well as the following timc-indcpcndcnt rela- 
tions between F, the power law index a and the equation 
of state parameter [l^, [l^l i 



wb — 2(F — 1) awB — 
1 + 2(r - 1) ~ a + 2 



where F = 1 -I- i > 1 from Eq. ([6]) for the IPL potential, 
and Wb is the equation of state of the fluid component 
dominating the background. So, during the era of radia- 
tion domination, with wb = \, 
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and during the era of matter domination, with, wb = 
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We also note here that, as in the case of all tracker poten- 
tials, the tracker solution for the IPL model is approached 
differently for different initial conditions. For example, 
in what is referred to as the "overshoot" case, p^j begins 
from a value greater than the tracker solution value. As- 
suming that (j) is released from rest, the dynamics of the 
scalar field start with an early kinetic phase (0^ » V) 



in which ui ^ 1 so that 



(12) 



c a^^ (from Eq. and V 

decreases very rapidly as <f) runs downhill. Since the ki- 
netic energy is too large for cf) to join the tracker solution 
as (f) rolls further down the potential, (j) will overshoot the 
tracker solution. The field will then freeze (as will V and 
■^) as rushes towards —1. Finally, when (f) rejoins the 
tracker solution, 4> will run downhill again and will 
increase from -1, briefly oscillate, and then settle into the 
tracker value p^ . 

In the "undershoot" case, begins from a value 
much smaller than the tracker solution value, and 4> is 
once again released from rest. This corresponds to the 
kinetic energy density being very small and </>, V , and 
being approximately constant or "frozen" as the universe 
evolves. Then, as in the "overshoot" case, reaches 
close to —1, p0 ~ const., and ps is decreasing. The 
value of W0 then increases from —1 as once again runs 
downhill. After a few oscillations, will then rejoin the 
tracker solution until p^ becomes the dominant compo- 
nent in the universe. 

Figures [H O and [J] depict the evolution of for the 
IPL model during these various regimes. With little sen- 
sitivity to the exact value of Vq, a will determine the 
amplitude of the w curve and determine the value of 
Wo = w{z = 0) > — 1 as long as « Mp. For given 
values of a, 4>i determines when the scalar field joins the 
tracker solution and how long it follows the tracking so- 
lution (Fig. [3]). As is pointed out in [2^, we also find 
that for the smaller values of a that we focus on in this 
work (e.g., a < 1), the smaller a is, the later the tracker 
is reached for a given initial value of (j) (Fig. 2]). With re- 
gards to Vb, we find that while increasing (decreasing) the 
value of Vq leads to corresponding increases (decreases) 
in ujde ~ y'^^ ^=0 (where h = ^), as expected, it 
leads to very small (essentially negligible) decreases (in- 
creases) in the value of wq and essentially no change in 
the tracking solutions or tracking behavior. When the 
scalar field has tracking solutions, different values of (f>i 
lead to similar values of, for example, —0.9 > wq > — 1, 
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with w — > —1 and — > as a — > oo. There wih be 
essentially no dependence of 0/ << Mp on either the 
present dark energy equation of state or the present con- 
tribution of dark energy to the total energy density of 
the universe (as illustrated in Fig. [3]). 
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FIG. 2: An Illustration of how the evolution and tracking 
behavior of w as a function of scale factor a is affected by 
different values of a, Vq, and (pj. The a scale is logarithmic 
here in order to show behavior on all time scales. 
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FIG. 3: Examples of how the evolution and tracking behavior 
of w as a function of scale factor a is affected by different 
values of (f)i for given values of Vo and a. For all curves, Vb = 
0.38 and a = O.I. These examples illustrate how different 
values of 0/ lead to the same values of the equation of state 
parameter today. The a scale is logarithmic here in order to 
show behavior on all time scales. 
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FIG. 4: This figure depicts the evolution and tracking behav- 
ior of ui as a function of scale factor a for different values of 
a for given values of Vo = 0.38 and cf)i = 10""^". As long 
as « Mp, a will determine wo and the amplitude of 
the w{a) curves. In addition, the smaller a is, the later the 
tracker is reached for a given 0/. The a scale is logarithmic 
here in order to show behavior on all time scales. 



solutions do not strictly exist and thus tracking behav- 
ior does not strictly occur for, roughly, all 0/ > 10~^ 
when a < 1 and 0.25 ^ Vq < 0.45. Moreover, for some 
instances in which —1.5 < logiQ((j)i) < —0.3, w « — 1 
initially but then increases towards —1 > wq > —0.9, for 
exarnple, as for the case of "thawing" models and behav- 
iors [23, US- Examples of this nontracking "thawing" - 
like behavior of the equation of state for 0/ = 10""'^ for 
a = 0.2 and 0.1 can also be seen (dashed-dotted curves) 
in Figs. [2] and [3l Nontracking initial conditions for the 
IPL model as well as possible connections between the 
quintessence field and the inflation field (the inflaton), 
which is beyond the scope of this work, are discussed in 
some detail in [2^ and references therein. Like [25|, and 
as we discuss further in Section IIVBI we also find that 
for cases where (/)/ — > Mp and the field has not joined the 
tracker by the present epoch, the range of acceptable val- 
ues of a increases significantly as w increases. For values 
of logio{(j)i) roughly between —5 and —l,w^ leaves its 
tracking phase with matter and enters a transient phase 
(see Fig. [3]) before exhibiting "thawing" behavior for 
logio{4>i) > -1.5 



D. The transition from tracking to acceleration 



C. The non-tracking case 

It is possible to find non-tracking cosmological solu- 
tions for IPL quintessence. If </>/ ^ Mp, then will 
follow the tracker solution for only a very brief period of 
time or not exhibit tracking behavior at all. In our com- 
putational algorithms, for example, we find that tracking 



For most of this work, we focus on cosmological solu- 
tions that exhibit tracking at early times. Out of respect 
for big bang nucleosynthesis [2Q] and other standard con- 
siderations there must be an early epoch of radiation 
domination where << pr and redshifts as P, [27| 

4a 

P0 oc a 2+0. (15) 
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It is possible in this case to find an exact solution to the 
Klein Gordon equation for which 

4)(xa^, (16) 

and it can be shown that this solution is an attractor [l[ . 
During matter domination, the attractor is also charac- 
terized by the scalar field evolving as 

(j)(xa^, (17) 

corresponding to energy density evolving as 

/90cxa~^. (18) 

As long as ^ << 1, these expressions provide a 
very good approximation to the behavior of the IPL 
quintessence field [H, In other words, the track- 

ing regime itself is strictly valid only when the expansion 
of the universe is dominated by matter. Then, at later 
times, when starts to make a significant contribution 
to the cosmic expansion rate, the value of Wrj, in Eq. (|12|) 
starts to diverge from its tracker value, as do 4>{a) and 
P4,{a), such that the scalar field mimics the behavior of 
a cosmological constant today (with w « —1), consistent 
with current observations. So, we can see that in the 
attractor solution decreases less quickly than pm and pr, 
which allows us to realize the following behavior: Deep 
in the era of radiation domination, p^ is small enough 
to satisfy constraints from standard models for big bang 
nucleosynthesis and the formation of the light elements, 
but p^ does eventually become large enough today (with 
w ^ — 1) so that the universe undergoes accelerated ex- 
pansion and acts as if it has a cosmological constant, but 
one that slowly varies with time and position [29| . 



argued that such small values of a lead to smaller basins 
of attraction and thus some degree of fine-tuning and 
dependence on initial conditions for the IPL model. We 
have observed, however, that for the realistic cosmologies 
that we consider for this work there remains a substantial 
basin of attraction: We can vary the initial conditions 
over a very large range of values with the end results for 
fl^fi, for example, still being physically acceptable [43 |. 

Other authors (e.g.. [H, [13, SI]) have also explored 
a variety of issues associated with tracking properties 
and solutions for this model. They considered theoretical 
constraints relating to, for example, equipartition initial 
conditions between quintessence and the remaining fiuid 
components which argue for larger values of a [isl . [25j . 
However, in our work we have focused for the most part 
on realistic families of cosmological solutions that are 
broadly consistent with observational constraints (i.e., 
a < 1) and which also include IPL tracking properties 
and behaviors that give the model its conceptual appeal. 
We also note that [l3| have found that while a is tightly 
constrained, IPL models with 0.25 < flm ^ 0.4 remain 
viable. 

The real appeal of IPL models from our point of view is 
that they offer an interesting class of non-A cosmologies 
with some degree of theoretical motivation. Thorough 
discussions of the basin of attraction (as well as the still 
outstanding cosmological constant problem) are key to a 
fundamental understanding of the ultimate importance 
one might give to the IPL model. We regard such dis- 
cussions as too poorly developed at this point to give 
them much weight in the very phenomenological analysis 
in this paper. For our purposes, it is good enough that a 
large range if initial conditions can converge to a common 
solution thereby avoiding to a substantial degree the fine 
tuning of initial values of — and [25j . 



E. Current constraints 

From an observational standpoint, if we require wq 
to be roughly consistent with current observational con- 
straints, say, for example, — 1 < w < —0.8, [sollsil. [s^. [ssj 
then the power law index a must be roughly in the range 
^ ci ^ 0-5, yielding a shallow potential shape. The 
quintessence equation of state in the current epoch aban- 
dons the tracking regime because the dark energy is now 
the dominant component. However, the shallow poten- 
tial shape makes wq not far from the tracking one in 
Eqn. nil), differing typically at the 10% level [2i|,[3l|. 

Various combinations of data (including CMB and SNe 
la observations) have been used to constrain the slope of 
the IPLpotential, finding a < 1 - 2 (e.g., [Hill 111,113, 
m, Sliyl), so that flatter potentials seem to be favored 
by the data. Recently, for example, [1^ have found 0.7 < 
a < 0.8 in an MCMC analysis of the IPL potential when 
assuming that the energy scale of the potential is that 
of a cosmological constant (i.e., Vq ~ A^ « 
and for when = varies in the range 0.1-0.9 and 
h = 0.70. A number of authors (e.g., [H, |4l|, \^) have 



III. PARAMETERIZATION OF THE INVERSE 
POWER LAW MODEL 

As a general rule, MCMC analysis requires a care- 
ful choice of the model parameters to be varied. Poor 
parameter choices and degeneracies between parameters 
can slow the rate of convergence and mixing of the 
Markov chain, reducing the overall efficiency by which 
the Markov chain explores a parameter space. For the 
IPL potential, V = Vb(^^)", the obvious choice of po- 
tential parameters to be varied is (/)/, a, and Vq. When we 
carried out our MCMC analysis of data forecast by the 
DETF to constrain the IPL quintessence model around a 
fiducial ACDM model, we chose our fiducial value for Vq 
(in units of h"^) to be 0.38, which is the value of the dark 
energy density today for a cosmological constant. We 
chose to make Vq a model parameter in our MCMC anal- 
ysis rather than keeping it fixed because other choices of 
Vq could provide equivalent cosmological solutions, and 
we were also interested in ascertaining how the MCMC 
exploration of the parameter space and its ability to con- 
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strain the other parameters would be affected by varying 
Vq as well. 

We have not found a need to reparameterize the IPL 
parameters to the extent that has been done, for example, 
in [1,0] for the Albrecht-Skordis or Exponential potential 
quintessence models. We did, however, find it necessary 
to place bounds on some of the potential parameters in 
order to prevent the MCMC from infinitely stepping into 
divergent directions of parameter space and thus never 
converging to a stationary probability distribution. An- 
other reason we placed bounds on the potential param- 
eters was to prevent the MCMC from spending possibly 
large amounts of computer time exploring uninteresting 
regions of parameter space that may be completely incon- 
sistent with observational and theoretical constraints. 

We placed a lower bound of on a, as a > is required 
for the pure IPL model that we consider . Given that 
the DETF data used in the first part of our MCMC anal- 
ysis is modeled around a cosmological constant, the most 
probable values of a will be those in which a approaches 
zero. From Eq. ([7]) we see that as 0/ Mp any value of 
a will lead to the same value of the potential V{4)) for a 
given Vq. However, since a largely controls the shape of 
the potential (as well as the amplitude of w{a)) and thus 
the evolution of the dark energy density and w^^, we find 
that the simulated data sets place sufficient constraints 
on a to prevent the MCMC from infinitely stepping into 
divergent directions in the a — (pi and Vq — 0/ parameter 
spaces even when (/>/ —>■ Mp. This renders a stringent 
upper bound on a unnecessary. 

We can also see from Eq. ([7]) that 0/ can take on any 
value and lead to solutions indistinguishable from a cos- 
mological constant as a — > 0. This degeneracy leads to 
a divergent direction in the a — space, where 4>i can 
be arbitrarily large or small. Also, the simulated data 
sets do not constrain (pi nearly as tightly as a due to 
the fact (previously discussed in the context of attractor 
solutions) that a broad range of (pj values can lead to the 
same 0o f^nd wq and thus have little effect on the evolu- 
tion of the dark energy density. Because of this effect, it 
is necessary to choose some cut-offs on (pj so that these 
infinite directions are bounded. 

As discussed in Section FlIB I we have parameterized our 
potential in a way that gives cosmologically realistic so- 
lutions where V{p) approaches the value of the dark en- 
ergy density today when p « Mp. With this in mind, 
we impose an upper bound of Mp on (pi which helps 
avoid solutions with uninterestingly low values of as 
well as solutions that are dominated by transients. We 
also note here that, given that the main thrust of our 
work involved an MCMC analysis of the regions of pa- 
rameter space associated with tracking, we have selected 
or filtered out non-tracking parameter values in the al- 
gorithms used to generate likelihood contours from the 
MCMC chains by implementing in our algorithms the cri- 
teria for tracking solutions (as discussed in Section fll Al 
and III Bp and, specifically, the "equation of motion" dis- 
cussed in [13]. Thus, all of the error contours displayed 



and discussed in sections IIVBI and IIV CI correspond to 
portions of the parameter space associated with tracking 
(i.e., parameter values corresponding to attractor solu- 
tions of the Klein-Gordon equation). Incidentally, we 
have found that for a typical Stage 2 MCMC chain gen- 
erated from a KCDM model, for example, about 90% of 
points stepped to in the chain correspond to parameters 
with tracker solutions, whereas the other 10% correspond 
to non-tracking (transient and thawing) parameters. 

Regarding a lower bound on </>/, we recall from Section 
IIIBI that we must have (pi « Mp so that the present 
field value, p^ (of order Mp), is reached from a very broad 
range of initial conditions. This insures that the tracking 
properties and solutions that make this model appealing 
are still included and valid within the parameter space 
explored in our MCMC analysis. If the lower bound on 
(pi is too large (pi may reach the tracking phase only 
at very late times or only by the present time (or not at 
all) , leading to a small basin of attraction and fine-tuning 
problems. We find that placing a lower bound of 0/ = 
10^^" in our MCMC analysis gives reasonable results by 
ensuring that on the one hand the tracking solutions and 
properties are included in the parameter space explored 
by the MCMC (i.e., there is a larger basin of attraction 
and (po fa Mp) but on the other hand, an appropriate cut- 
off or bound has been placed on a divergent direction in 
the a — (pi space that may not otherwise be constrained 
by the data (and thus possibly preventing the MCMC 
chains from coming to equilibrium). 

The above lower bound is not well suited for examining 
the finer details of nontracking transient and "thawing" 
regions of parameter space (where (pi Mp). In chains 
with a lower bound of 10"^" or smaller on (pi, the part 
(« 10%) of the chain that shows nontracking and thawing 
IPL solutions is not sufficiently well populated to show 
the full structure of the probability distribution. In or- 
der to allow the MCMC to step more frequently in these 
parameter space regions and so better converge (as dis- 
cussed in [1|) on a well-resolved probability distribution 
for the nontracking and thawing regions of the parame- 
ter space, we have also carried out an MCMC analysis 
with a lower bound of —3 placed on logio{(pi) (see Section 
HVDll . 



IV. MCMC RESULTS AND ANALYSIS 

A. General approach 

Following the approach taken by the DETF, we gen- 
erated "data models" or simulated data sets for future 
SNe la, baryon acoustic oscillation (BAO), weak gravita- 
tional lensing (WL), and CMB (PLANCK) observations. 
These considerations of DE projects follow developments 
in "stages" : Stage 2 represents ongoing projects that are 
relevant to dark energy; Stage 3 consists of medium-cost, 
near-term, currently proposed projects (such as BAO, 
SNe la, and WL surveys with 4-meter class telescopes 
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using photometric redshifts); Stage 4 consists of a Joint 
Dark Energy (Space) Mission (JDEM), Square Kilometer 
Array (SKA), and/or Large Survey Telescope (LST) 0. 
"Optimistic" and "pessimistic" versions of the same data 
models give different estimates of systematic errors. Ad- 
ditional information on the specific DETF data models 
is given in Appendix A of [3| and the technical appendix 
of the DETF report [1]. We excluded the DETF galaxy 
cluster data models in our work because the extension 
of the DETF calculations to our analysis is not straight- 
forward, especially in regards to estimates of systematic 
errors flSSIil. 

We have generated two sets of data models. One type 
is generated around a cosmology with a cosmological con- 
stant, consistent with DETF Stage 2, 3, and 4 SNe la, 
WL, BAO, and CMB data models. The other set of 
data models is built around an IPL fiducial model which 
was chosen to be consistent with simulated Stage 2 data 
based on a cosmological constant cosmology. We then use 
an MCMC algorithm to map the likelihood around each 
fiducial model [KCDM and IPL) via a Markov chain 
of points in parameter space, starting with the fiducial 
model and moving to a succession of random points in 
space using a Metropolis-Hastings stepping algorithm. 
The technical details of our MCMC algorithm are pre- 
sented in Appendix B of Q and references therein. In 
this way we can, for example, analyze the parameter 
space of IPL quintessence in the light of DETF data 
models and evaluate the likelihood function of the pa- 
rameters of our model. Once the Markov chains of our 
models in parameter space have been computed we can 
extract likelihood contours from the distribution of mod- 
els and display them as projected 2-D likelihood contour 
plots. This can then give us a picture of the shape of the 
likelihood region of all the parameters in our models in 
the whole multidimensional parameter space if we were 
to plot likelihood contours for each pair of parameters 
in the parameter space. In all plots in this paper, we 
show 68.27% (Icr), 95.44% (2cr), and 99.73% (3cr) confi- 
dence contours, which consist of points where the likeli- 
hood equals e~^2~, e~^^ , and e 3- of the maximum 
value of the likelihood, respectively. We have constructed 
these plots by marginalizing over all of the cosmological 
parameters, uSm, ijJk, ^b, h, S^, ns, n'^, (as defined by 
the DETF), and the various nuisance and/or photomet- 
ric redshift parameters, which take into account uncer- 
tainties and errors in the simulated data. The nuisance 
and photometric redshiftparameters are described and 
explained in detail in [j, Q . 



B. Cosmological Constant Fiducial Model 

In this section we present the results of our MCMC 
analysis for the combined simulated data sets generated 
around a ACDM cosmology. We list the values of the free 
parameters for our ACDM fiducial model (with energy 
density and Vq in units of and ipi in reduced Planck 



units) in Table [H (The IPL parameters given generate a 
cosmological constant.) 

TABLE I: Fiducial Parameter Values (energy densities in 
units of h^) for ACDM model. 





0.3796 


OJm 


0.146 


(^k 


0.0 


LOB 


0.024 


LUr 


4.16 X 10 


Us 


1.0 


n's 


0.00001 


h 


0.87 


h 


0.72 


a 


0.0 






Vb 


0.38 



We note that h^ia = 1) = ojm -\- Wr + uJk + (jJoe, with 
recent observations providing a prior constraint of ft- = 
0.72 ± 0.008 dil. Also, w^, the radiation energy density, 
is not a free parameter for our calculations but is fixed 
by the CMB temperature (and the standard assumption 
of three massless neutrinos) 

Stage 2 combines SNe la, WL, and CMB data mod- 
els but does not include BAO data models. Stages 3 
and 4 additionally include the BAO data models as well. 
As discussed in Section flV Al we project our probability 
distributions into 2-D spaces given by pairs of the IPL 
parameters (i.e., the Vb — a, Vb — 0/, and <f>i — a planes). 

The likelihood contours in the ~ a plane, with all 
non-tracking parameter values {logiQ{(f)) < —6) excluded, 
for Stage 2 and the optimistic versions of Stage 3 pho- 
tometric, Stage 4 Space, and Stage 4 Ground LST com- 
bined data are shown in Fig. \5\ In all cases the error 
contours show the expected trend of the IPL potential 
to approach a cosmological constant as a —> (and also 
corresponding to where the slope of the potential goes 
to 0). The vertical axis where a = corresponds to A. 
Therefore, the value of ^(^o) = Vo on the vertical axis 
represents A or the dark energy density lude for a = 0. 
However, along the lines of the discussion in Section III 
of for the Albrecht-Skordis model and as discussed in 
this paper in Section FlIBl we must also keep in mind that 
the parameter Vb does not have a significant effect on the 
equation of state of dark energy. Moreover, for a > 0, Vb 
is no longer identical to lude,o = ^de{z = 0). 

For small values of a, there is a spread in Vq in the 
Vb — a space. Since these values of a are consistent with 
A or a non-evolving dark energy, the spread in Vq is es- 
sentially a measure of how well the experiments are mea- 
suring ujde.o- The spread or uncertainty in Vb for all 
a is also a result of uncertainties on measurements of 
ri„i_o- Larger values of a correspond to larger values of 
w {w > —1) and thus values of w that deviate more and 
more from the equation of state for A as a increases, 
possibly up to values of a that correspond to detectable 
differences from A. The smallest values of Vb correspond 
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FIG. 5: Vo-a lcr(68.27%), 2cr(95.44%), and 3cr(99.73%) con- 
fidence regions for DETF "optimistic" combined KCDM data 
models. 



reduction in the Vq direction with increasing stage num- 
ber, indicating the improving constraints that the data 
places on the dark energy contribution to the total energy 
density of the universe today. 
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to the largest values of a, which in turn correspond to the 
largest values of the equation of state (and hence those 
values of w deviating the most from what we expect for 
a cosmological constant). As the value of a increases, 
we see that the likelihood contour in Fig. [5] has an over- 
all downward curved shape. This is due to the fact that 
the slope of the potential becomes steeper for increasing 
values of a, which leads to greater evolution of the dark 
energy density and larger values of wq that deviate more 
and more from —1. The reduction in the Vb direction re- 
flects improving constraints with increasing stage number 
on the dark energy density. As a specific quantitative ex- 
ample of this, we see from the Stage 4 error contours in 
Fig. [5] that the extrema of the range of Vb values deviates 
from the fiducial value by less than 20% when a « 0.1 
and less than 5% when a w 0. The shrinking in the a 
direction corresponds to increasing constraints on devia- 
tions from a cosmological constant. 

Fig. [HI depicts likelihood contours in Vb — logiQ{(j>i) 
space, where, again, all non-tracking transient and thaw- 
ing parameter values have been removed. As noted in 
Section Hill we imposed 10"^" < (j)i/Mp < 1. Since a 
large range of initial values of the scalar field (</)/ < Mp) 
are generally washed out by the tracking behavior, we 
can see from the contours that there is very little depen- 
dence of the dark energy density today on 0/ << Mp. 
Once again the spread in Vb values is essentially a mea- 
sure of how well the experiments are measuring the dark 
energy density at the present time. The error contours 
also show a slight trend toward an increasing range of 
acceptable values of (pi which possess attractor solutions 
as Vb decreases, which is associated with greater a values 
and thus greater dark energy evolution. The sections of 
the overall parameter space depicted in these figures also 
tend to disfavor larger values of a, or, equivalently, dis- 
favor larger departures from a cosmological constant and 
thus more dark energy evolution. We once again note a 



FIG. 6: Vo - logio{<j)i) lcr(68.27%), 2cr(95.44%) and 
3(t(99.73%) confidence regions for DETF "optimistic" com- 
bined KCDM data models. 
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FIG. 7: logw(4>i) ~ a lcr(68.27%), 2cr(95.44%) and 
3(t(99.73%) confidence regions for DETF "optimistic" com- 
bined KCDM data models. 

The likeHhood contours in the log{4>i) — a (Fig. [7]) 
space are clearly seen to shrink in the a direction with 
increasing stage number, once again showing improving 
constraints on the amount of dark energy evolution and 
on deviations from a cosmological constant from Stage 2 
to Stage 3 and from Stage 3 to Stage 4. This corresponds 
to a greater disfavoring of larger values of a with succes- 
sive stages of data. We also see in Fig. [7] a very slight 
trend toward an increasing range of acceptable values of 
(j)i possessing attractor solutions as a increases. This cor- 
responds to the trend of a larger range and upper limit 
for (/)/ having attractor solutions for smaller values of Vb 
discussed in regards to Fig. [HI 

The trend in Fig. [7] is related to the fact that the 
largest values of of (pi from which the attractor is joined 
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before the present time occur on the flatter portions of 
the potential where is closer to —1 in recent times 
and today and the curvature and slope of the potential is 
smaller. Since a controls the steepness of the potential, 
changes in a have less of an effect on the flatter parts 
of the potential where </>/ is larger and V{(j)i), V^((/)/)', 
and V{4>i)" are smaller (as can also be seen in Eq. 
Eq. (HH])). So, when the scalar fleld tracks the background 
evolution on flatter portions of the potential, we expect 
a slight increase in the range of acceptable (tracking) (pj 
values as a increases. 

Overall, as found by the DETF, successive stages of 
data do better at constraining the evolution of dark en- 
ergy. As can be seen in the likelihood contours above and 
as was also found for the case of the Albrecht-Skordis 
model 0, the IPL potential parameters appear to be 
somewhat better constrained by the DETF Stage 4 LST 
ground data models than by the DETF Stage 4 space 
data models. This reflects the fact that ground and space 
data are sensitive to slightly different features of the dark 
energy evolution. 



C. Inverse Power Law Fiducial Model 

We next evaluate the power of future experiments by 
assuming that the dark energy in the universe can actu- 
ally be described by the inverse power law model rather 
than a KCDM fiducial model. For our fiducial IPL 
model, we use a = 0.14, 0/ = 10~^^, and Vb = 0.31. 
The remaining parameters of the IPL fiducial model are 
the same as those used in the fiducial KCDM model. 
Our IPL model fiducial values (given in Table |TT] with 
energy densities and Vq in units of h? and 4>i in reduced 
Planck units) were chosen, excluding consideration of the 
"thawing" or outlying regions of the parameter space, to 
lie near the boundary of (or just beyond) 1 a detection 
or within the 95.44%(2(t) confidence region in the Vb — a 
and log{(j)j) — a spaces (Fig. [5] and Fig. [7]) for Stage 2 
AC DAI data, but excluded by more than Scr in the Stage 
4 optimistic ground and space data so as to be strongly 
ruled out by Stage 4 KCDM data. 

TABLE II; Fiducial Parameter Values (energy densities in 
units of h^) for Inverse Power Law model. 
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We also ensured that this fiducial model had initial 
conditions and had an equation of state such that the 
attractor is joined before the present time. The equation 
of state parameter as a function of scale factor a for all 
time scales (the a scale is logarithmic) for our fiducial 
model is similar to the dashed curves in Figs. [2] and 
m We also depict the potential of the fiducial model in 
the top panel of Fig. [5] along with the corresponding 
equation of state evolution as a function of redshift in 
the bottom panel. The fiducial model corresponds to 
the point wq = —0.955, which deviates from u;(z) = — 1 
by only about 4.5%. We have chosen our fiducial model 
to thus be marginally consistent with the ACZJM-based 
data but demonstrating enough dark energy evolution to 
be different enough from A to be resolved by Stage 4 
experiments. In this way we are able to illustrate the 
power of Stage 4 data models and their ability to rule 
out the A model. 
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FIG. 8: The potential of the IPL fiducial model (a = 0.14, 
(l>i = 10"^^ Vo = 0.31) (top panel,dashed curve). The cor- 
responding equation of state evolution w{z) for a potentially 
observable range of redshift values is shown in the bottom 
panel. The solid curve overlaying the potential in the top 
panel shows the evolution of the IPL fiducial model scalar 
field for the range of z values (from 2 = 5 to the present 
time) depicted for w{z) in the bottom panel. 

Duplicating our MCMC analysis methods for the IPL 
fiducial model, we again marginalized over all but two 
pairs of the parameters a, ipi, and Vq for the purposes 
of generating 2-D likelihood regions for the IPL dark en- 
ergy parameters. Fig. [9] shows the results of our MCMC 
analysis and calculations for Stage 2, Stage 3 Photo- 
optimistic, Stage 4 LST Optimistic, and Stage 4 Space 
Optimistic data models in the Vq — a parameter space. 
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We can see from the a = axis, corresponding to a 
cosmological constant, that the AC DM model (i.e., a 
non-evolving scalar field) is still allowed at Stage 2 (at 
the 2(7 (95.44%) confidence level but not quite at the 
1(7 (68.27%) confidence level) but becomes less favored 
by subsequent stages of data models. At Stage 3 the 
ACDM model lies outside of the 2a contour, and by 
Stage 4 it is ruled out by well over 3cr. For Stage 2 and 
subsequent stages the range of a values covered by the 
contours is significantly greater than for the ACDM case 
since dark energy solutions with more evolution are fa- 
vored more here. The greater dark energy evolution for 
this case also leads to the slightly more significant down- 
ward trend in the shape of the contours than is seen in 
the AC DAI confidence contours. 
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FIG. 9: Vo-a la (68.27%), 2a (95.44%) and 3a (99.73%) like- 
lihood contours for DETF optimistic combined data sets gen- 
erated from a selected IPL background cosmological model. 



The described increase in constraining power for higher 
quality data models is similar to the ACDM results in 
Section flV Bl for the ACDM model. However, as previ- 
ously indicated, the range of a values has significantly 
increased within the la, 2a, and 3a contours, allowing 
for an increased range of evolving dark energy solutions. 
By the Stage 4 combined data sets, we can clearly differ- 
entiate between our selected IPL fiducial model and the 
ACDM model by well over 3a. This increased constrain- 
ing power is again consistent with the {ACDM) DETF 
results for Stage 4 experiments. Hence, the results of our 
MCMC analysis, as seen in Fig. [5] (as well as Fig. [11] be- 
low), show that, for a universe described by this specific 
IPL fiducial model, the Stage 4 experiments will rule out 
a cosmological constant by well over 3a. 

Figure \W\ shows likelihood contours in Vq — log{4)i) 
space. As for the case of the ACDM model, there is 
again very little dependence of dark energy density to- 
day on 07 when 0/ << Mp. Once again, the spread in 
Vo values is essentially a measure of how well the exper- 
iments are measuring the present dark energy density as 
given by the chosen IPL fiducial model. The trend to- 



ward an increasing range of and acceptable upper limit 
to values of </>/ possessing attractor solutions for smaller 
Vo, as noted in reference to Fig. [SI is sHghtly more pro- 
nounced here due to the larger range of acceptable values 
of a and greater DE evolution for the IPL model. 
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FIG. 10: Vo 



la (68.27%), 2cr (95.44%) and 3cr 



(99.73%) likelihood contours for DETF optimistic combined 
data sets generated from a selected IPL background cosmo- 
logical model. 



The likelihood contours in the log{(f)i) — a (Fig. [TT|) 
plots are clearly seen to shrink in the a direction with in- 
creasing stage number, but the overall range of a values 
stepped to by the MCMC chain and thus included within 
the likelihood contours is significantly larger than for the 
ACDM model data sets, again indicating that dark en- 
ergy solutions with more evolution are disfavored less for 
this IPL fiducial model than for the ACDM model. The 
log{<j)i) — a contours also show (like the Vo ~ log{(j)i) con- 
tours in Fig. [To]) that the ACDM model is still allowed 
at Stage 2 (but lies just outside of the la contour) and 
at Stage 3 Photo-optimistic (lying outside of the 2a con- 
tour here) but is ruled out by well over 3f7 by Stage 4, 
again becoming less favored by subsequent stages of data 
sets. We also see a slightly more pronounced trend of an 
increasing range of acceptable values of 0/ possessing at- 
tractor solutions as a increases. Again, this corresponds 
to the increasing range of acceptable 4>i values possess- 
ing attractor solutions for smaller Vq in Fig. [TO] and the 
fact that the part of the IPL potential where the largest 
values of 0/ that lead to attractor solutions that are still 
acceptable is steeper (larger a) than for the ACDM case 
(Eqn. ((S])-Eqn. ([TO])). 

As in the case of the ACDM data sets discussed in Sec- 
tion IIVBI and was found in the Albrecht-Skordis model 
[(|] , we find once again that Stage 4 ground data (this time 
based on our fiducial IPL model) slightly more strongly 
constrains the parameters a and Vq than does the Stage 
4 space data. This is opposite of what has been found 
with other scalar field models [1, Q . 
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FIG. 11: log{(l>i) - a la (68.27%), 2a (95.44%) and 3a 
(99.73%) likelihood contours for DETF optimistic combined 
data sets generated from a selected IPL background cosmo- 
logical model. 



D. Non- Tracking Parameter Space Regions 

Though the main focus of our work has involved an 
analysis of the tracking regions of the parameter space 
of the IPL model, here we discuss briefly the results of 
our MCMC analysis of the non-tracking regions, i.e., ini- 
tial values of the scalar field from which the attractor 
is not joined before or by the present time. In this case, 
our motivation is simply to explore an interesting-looking 
class of dark energy behaviors that have already been 
considered elsewhere in the literature (e.g., [1^). We ac- 
knowledge that to the extent that the tracking behavior 
is a key reason to consider the IPL model, the solutions 
considered in this section do not benefit from the same 
degree of motivation. 

As indicated previously, in order to allow the MCMC 
to step more frequently in non-tracking portions of the 
parameter space and so bring out greater detail in the 
thawing and some of the transient portions, we have also 
generated MCMC chains with a lower bound of —3 placed 
on logi{){4>i). These outlying regions of parameter space 
associated with the thawing equation-of-state behavior 
(again corresponding to 4>i Mp and increasing w and 
present-day dark energy density values in recent times) 
can clearly be seen in the Stage 2 and Stage 3 error con- 
tours for our IPL fiducial model depicted in Fig. [T^] at 
a > 0.5, where the contours turn or "flare" upward and 
become more "patchy" . This un-smooth and flared ap- 
pearance of the 2 and 3tT contours correspond to the 
largest values of 0/, where the equation of state w{a) 
increases or does not turn down as steeply near scale fac- 
tors of unity and so is exhibiting thawing-like behavior, 
and, therefore, the acceptable range of a values signifi- 
cantly increases as w increases. 

These portions of the likelihood contours correspond 
to outlier points lying relatively far outside the main dis- 
tribution of parameter points stepped to by the MCMC 



chain. In these regions of parameter space the scalar field 
starts to evolve on the fiatter portions of the IPL poten- 
tial where V{(l)i) is small. We see that there is a greater 
spread in Vq values, and, thus, Vq is less constrained by 
the data here. This is related to the fact, again, that the 
corresponding equation of state values for (f>j ^ Mp don't 
turn down as steeply near scale factors of unity (or even 
increase towards values greater than —1) compared to w 
values corresponding to (pi < Mp. We can see that the 
area of the these outlying likelihood contours shrinks and 
tightens from Stage 2 to Stage 3 and again from Stage 3 
to Stage 4. The reduction in the Vq direction again shows 
improving constraints with increasing stage number that 
the data places on these outlying transient and thawing 
regions. We also observe an apparent illustration here of 
the ability of the Stage 4 ground-based simulated data 
sets to better constrain the thawing behavior than the 
Stage 4 space-based data. This appears to be consistent 
with the results obtained in the MCMC analysis for the 
tracking regions of parameter space. Moreover, the Stage 
4 space-based data rules out a significant portion of the 
thawing region of the parameter space, while the Stage 
4 Ground LST Optimistic data sets appear to rule out 
nearly all of the thawing parameter values. 
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FIG. 12: Vo-ala (68.27%), 2a (95.44%) and 3a (99.73%) 
likelihood contours for DETF optimistic combined data sets 
generated from a selected IPL background cosmological model 
for the case of a cut-off of logio{4ii) = —3 placed on the 
MCMC algorithm. This effectively gives an enlarged and 
more detailed view of non-tracking and "thawing" -like regions 
of the parameter space. 

Figure [T51 depicts log{(t)j) — a likelihood contours for 
our IPL fiducial model for non-tracking regions of param- 
eter space associated with transient and outlying thawing 
equation of state behavior. Given that our MCMC anal- 
ysis did not focus nearly as much on non-tracking regions 
of parameter space than the tracking regions, our chains 
may not have equilibrated for the non-tracking regions to 
the same extent as they have done for tracking regions. 
However, we believe that important trends can still be 
ascertained from this analysis. In the Stage 2 and Stage 
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3 likelihood contours, a significant increase in the a di- 
rection for the largest (/>/ values (</)/ Mp) can be seen. 
This corresponds to the fact, as discussed in Section HTB] 
that for the IPL model the range of acceptable a values 
is largest for the largest initial scalar field values from 
which the attractor is not joined by the present time. 
This is associated with the flatter part of the IPL poten- 
tial that is less sensitive to a, which controls the slope of 
the potential. 

As can be seen from Eq. ([9]), flatter parts of the poten- 
tial correspond to cases where V{(j)i) ~ V{(j)Q) is smaU 
and is large. Thus, even large values of a can be asso- 
ciated with flatter portions of the potential here. So, as 
long as V{(j)i) ~ V{4io) remains very small and the ratio 
of a and (/)/ does not become too large, a larger range 
of acceptable values of a, leading to similar cosmologies, 
will be allowed within the parameter space. Moreover, 
larger 0/ values combined with larger a can lead to sim- 
ilar w{a) > —1 with behavior close to that of A. Hence, 
for Stage 3 data and especially Stage 4 data, the MCMC 
will not step as much in this region (since AC DM mod- 
els and models with similar behavior are ruled out to a 
greater extent by Stage 3 and 4 data). This explains the 
greater constraints placed in the a direction for the very 
largest 0/ values for successive stages of data sets, and 
is consistent with the overall trend of the likelihood con- 
tours in the log{(j)i) — a space shrinking in the a direction 
with higher quality data. 

The extent to which larger a values (and thus signif- 
icant portions of the thawing regions of the parameter 
space) are constrained and even ruled out by the Stage 4 
data sets also reflects the degree to which evolving dark 
energy is constrained and disfavored by the higher quality 
data sets. Once again, and perhaps more dramatically il- 
lustrated here, we see that the ground-based Stage 4 data 
sets constrain the thawing regions of parameter space for 
the IPL model to a more significant extent than do the 
Stage 4 space-based data sets. 



V. DISCUSSION AND CONCLUSIONS 

We have presented our MCMC analysis of the inverse 
power law quintessence model using combined simulated 
data sets forecast by the DETF and representing future 
dark energy experiments. In doing so, we have analyzed 
the impact of DETF simulated data models in the con- 
text of the IPL model of dark energy and demonstrated 
the ability of these experiments to place significant con- 
straints on the parameters of a quintessence model. We 
have found that the effect of the DETF combined data 
models on the parameter space of IPL models is broadly 
consistent with the DETF findings. In particular, we 
have found a significant improvement in the constraining 
power of each successive stage of DETF simulated data 
sets. 

We have shown likelihood contours for choices of com- 
bined DETF data sets and found the increase in IPL dark 
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FIG. 13: log{(l)i) - a la (68.27%), 2cr (95.44%) and 3cr 
(99.73%) likelihood contours for DETF optimistic combined 
data sets generated from a selected IPL background cosmo- 
logical model for the case of a cut-off of log\o [4>i) = ~3 placed 
on the MCMC algorithm. This effectively gives an enlarged 
and more detailed view of non-tracking and "thawing" -like 
regions of the parameter space. 

energy parameter constraints with increasing data qual- 
ity to be consistent with the DETF results in the wo — Wa 
parameter space. For example, the relative constraints on 
the size of the Vq — 4>i parameter space between different 
simulated data sets lead to similar constraints computed 
by the DETF in the wq — Wa parameter space. A direct 
comparison with the DETF Figure of Merit was com- 
plicated by the fact that the IPL model depends on 3 
parameters (a, 0/, and Vb), while the DETF FoM was 
calculated based on the two-dimensional wq — Wa space. 
However, we found that the changes in the areas of pro- 
jected two-dimensional likelihood contours were consis- 
tent with the DETF results. Specifically, the DETF re- 
ported an FoM (defined as the inverse area inside the 95% 
likelihood contours in the wq — Wa plane) that showed a 
gain of at least a factor of 3 in going from Stage 2 to good 
combinations of Stage 3 data sets (and thus a factor of 
roughly 3 decrease in allowed parameter area when mov- 
ing from Stage 2 to good combinations of Stage 3 data), 
and a gain of at least a factor of 10 in going from Stage 
2 to good combinations of Stage 4 projects. We observed 
decreases by similar amounts in our projected 2-D likeli- 
hood contours for pairs of IPL parameters. 

In the course of this work we have also produced and 
examined similar 2-D likelihood plots of a much wider 
range of combined DETF simulated data sets, including 
data models with "pessimistic" estimates of systematic 
errors and data models representing single DE observ- 
ing techniques. We found our results in the IPL model 
parameter space to be consistent with the constraints re- 
ported by the DETF in the wq — Wa space across the 
complete range of data combinations and selections that 
we considered. 

We constructed our simulated data sets from two dif- 
ferent background cosmologies, one with a cosmological 
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constant and one with an IPL scalar field with specific 
parameter values. We found our results to be consistent 
with those of the DETF in both cases. We have sepa- 
rately analyzed cases constrained to having early track- 
ing behavior and other cases which focused on the non- 
tracking solutions. In each case we have placed bounds 
on some of the IPL potential parameters as necessary to 
prevent the MCMC from infinitely stepping in divergent 
directions of parameter space (and thus never converging 
to a stationary probability distribution) and to also en- 
able us to better examine and analyze details in enlarged 
regions of parameter space corresponding to non-tracking 
behavior. 

In order to demonstrate the power Stage 4 experiments 
will have for detecting the evolution of dark energy, we 
chose a specific background IPL scalar field model with 
parameter values of a = 0.14, 0/ — 10^^^, and Vq — 0.31 
that was consistent with Stage 2 data based on a cos- 
mological constant. This specific model corresponds to 
w{a = 1) = Wo = —0.95535, which deviates from w = —I 
by about 4.5%. One must look back to much earlier times 
(e.g., a < 0.2) and/or look to larger a parameter values 
in order to find more significant deviation from w = —1 
for this quintessence model (see Figs. [U [2l |3l HJ. We 
found that if the universe were in fact to be described 
by this fiducial IPL quintessence model, then good Stage 
4 experiments would rule out a ACDM model by bet- 
ter than 3(7, indicating that there is indeed a dynamical 
component to dark energy. For the IPL background cos- 
mology, we found that the ACDM model lies outside the 
la contour but within the 2a contour at Stage 2 and lies 
outside of the 2a likelihood contour by Stage 3. We also 
noted that the variable a was somewhat more strongly 
constrained by Stage 4 ground data sets than with Stage 
4 space data. This is consistent with the results reported 
by @ for a similar MCMC analysis carried out on the 
Albrecht-Skordis scalar field model, but is opposite of 
the behavior displayed by the Exponential and PNGB 
scalar field models as described in our other companion 
papers 0, Q . This effect is under current investigation 
and may lead to new insights into the complementarity 
of ground and space-based Stage 4 dark energy projects. 

We have found, as also discussed in H and demon- 
strated in our companion papers [1, H, 0], that widely 
varying families of functions w{a) for the IPL model are 
constrained by the DETF data sets in a similar way to the 
constraints found in the wq — Wa parameter space by the 
DETF. In particular, we have seen that the main IPL 
model potential parameter a is constrained by DETF 
data models in a comparable way to the constraints found 
in the wq — Wa formulation by the DETF, even though 



the Wo — Wa parameters describe very different functions 
w(a). We believe that this relates to the fact pointed out 
in [43] that high quality DETF data sets will be able to 
constrain many more properties of w{a) that are present 
in the wq — Wa parameterization alone and will thus be 
able to make good measurements of significantly more 
than two equation of state parameters. More specifi- 
cally, by considering the IPL family of w{a) functions and 
Wq — Wa family of functions in terms of an orthonormal 
basis of independently measure mode functions Wj (a) (as 
discussed in 0, S^), we are able to ensure that a wide 
variety of different w{a) functions will be constrained as 
well as the DETF wq — Wa parameters. In other words, 
the various quintessence models (discussed in this paper 
and in our companion papers) are just sampling differ- 
ent random combinations of the "well measured modes" 
discussed in [i^ and in each case lead to similar results. 
This also appears to reflect the fact that many more func- 
tions w(a) are measured than are contained in any of the 
quintessence model w{a) family of functions alone 
Consequently, modeling the impact of future dark en- 
ergy experiments using the two-parameter DETF scheme 
makes some sense in that it gives a good indicator of the 
impact of scalar field dark energy models with a similar 
number of parameters in the quintessence potential. 

One of the advantages of the techniques employed in 
this and the companion work 0, is that we can 

explicitly examine how simulated data sets representing 
future dark energy experiments can constrain actual the- 
oretically motivated quintessence models (in addition to 
abstract parameterizations such as the wq — Wa ansatz) 
in a significant way. As developed further in Y| this ap- 
proach helps us understand how future data has the ca- 
pability to reject some (or possibly even all) current dark 
energy models entirely. 
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